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1.0 INTRODUCTION 

1 .1 General 

This document presents the final report on the Shuttle/ 

TDRSS Modelling and Link Simulation Study performed for NASA 
Johnson Space Center under Contract NAS 9-14846 directed by 
William Teasdale and Sid Novosad. It represents a portion of 
the work accomplished during the period June 10, 1976 through 
April 20, 1979. 

The general objectives of the overall contract are the 
following: (1) to develop Shuttle/TDRSS link simulation models 
which can predict the performance sensitivity of these links to 
Shuttle transmitted signal distortions, (2) to exercise the simula- 
tion to predict the bit error rate performance as a function of the 
signal distortion parameters, (3) to present these data in a form 
suitable to support ESTL hardware tests and (4) to provide ESTL 
with a quick assessment of the basic differences between the ESTL 
TDRSS simulator hardware and the actual TDRSS system. This last 
item was covered in informal briefings to pertinent ESTL personnel 
and is not further addressed in this final report. 

In what follows an overall description of the contractual 
effort and a brief summary of the results is given. This is 
followed by backup material which includes simulation data and 
analyses from which our summary results and recommendations have 
been drawn. 

1.1.1 Final Report Contents 

This report addresses and documents LinCom's findings on the 
task statements detailed in the Statement of Work. All tasks 
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pertain to the performance prediction of Shuttle forward and 
return links through the Tracking and Data Relay Satellite 
System (TDRSS). 

Chapter 2 describes in detail the Shuttle/TDRSS S-band and 
Ku-band link simulation package, called LinCsim, developed for 
evaluation of link performance for specific Shuttle signal designs. 

The link models are described in detail and the transmitter distor- 
tion parameters, the so-called User Constraints are defined carefully. 
A complete list of the present values of all parameters used in the 
simulation is also included in the description. 

Chapter 3 contains the results of extensive simulation of the 
Shuttle/TDRSS S-band return link. The overall link degradation 
(excluding hardware degradations) relative to an ideal BPSK channel 
are shown for various sets of user constraint values and the perform- 
ance sensitivity to each individual user constraint is illustrated. 

The results can be used by NASA/JSC as a technical base for 
negotiating interface. parameters between the Shuttle and TDRSS. 

Chapters 4 and 5 take an indepth look at two specific problems. 
The effect of excessive Spacelab clock-jitter on the return link 
BER performance is analyzed in Chapter 4. This investigation 
accounts for the effects of the clock clean-up system (scrubber) 
in the Shuttle as well as for the coding and for the high level 
of data asymmetry expected on this link. .Chapter 5 considers the 
problem of subcarrier recovery for the K-band Shuttle return link 
signal when one of the channels contains unprocessed digital data 
from a detached payload. Two different implementations for the 
payload interrogator are investigated and their performance compared. 
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1 .2 Summary and Recommendations 

1.2.1 Shuttle S-band Return Link Performance Summary 

The performance degradation of the Shuttle/TDRSS return link 
due to Shuttle signal imperfections has been evaluated for the 
expected Shuttle distortion values based on current link budgets 
and TDRS/ground station hardware data. The results show a loss 
relative to ideal BPSK of 1.65 dB at the nominal bit error 
probability of 10"^. The TORS and ground station contribute 
.9 dB to this loss and the Shuttle the remaining .75 dB. This is 
illustrated in Fig. 1.1. 

The single biggest contributor to the Shuttle degradation 
is the phase noise which is specified as 10° rras in the frequency 
interval 0 to 270 Hz. We recommend to separate the phase noise 
into components lying inside and outside the carrier tracking loop 
bandwidth, respectively, since these components affect the bit 
error rate performance differently. Such a refinement of the 
specification could greatly improve the accuracy of the predicted 
performance results. Since our present prediction is based on 
worst-case assumptions the overall degradation would be reduced. 

1.2.2 Attached Spacelab Bit Jitter Performance Summary 

The functional model of the Spacelab return link is shown 
in Fig. 1.2. The data and/or clock transitions from the Spacelab 
into the Shuttle are subject to severe phase jitter. The clock 
regenerator (scrubber) removes some of this jitter on the reclocked 
and convolutionally encoded data. However, the coded data may 
exhibit up to 10% of data asymmeti^y. 
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Figure 1.2. Spacelab Return Link Block Diagram. 
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The effect of bit jitter in the scrubber was analyzed using 
measured pulse shapes. The results (Fig. 1.3] show that the error 
rate of the reclocked data is negligible for bit jitter values which 
can be tolerated on the Shuttle-to-ground link. This holds for the 
case where the mean sampling instant is in the center of the bit 
{0% bias in Fig. 1.3) as well as when it is offset by 12.5%, which 
is the maximum predicted value. It remains then to analyze the effect 
of the residual bit jitter and of the data asymmetry on the Shuttle- 
to-ground link performance. The analysis, based on an ideal clock 
in the ground station receiver, shows that the CNR loss due to 10% 
data asymmetry amounts to .52 dB for an otherwise perfect signal. 

The incremental • loss due to bit jitter is shown in Table 1.1. 

Based on spectral data provided by JSC personnel the clock 
phase jitter at the scrubber output was estimated to be less than 
2%. From Table 1.1 it is then clear that the associated CNR loss 
is well below .1 dB. 

1.2.3 Shuttle Bent-Pipe Mode Subcarrier Performance Summary 

A typical functional diagram of a detached payload return 
link is shown in Fig.. 1.4. The payload signal is received by the 
payload processor, converted to baseband, hard-limited and then 
modulated onto one quadrature component of the K-band return link 
subcarrier. The recovery of the subcarrier phase from such a signal 
was investigated and compared with the performance of a similar 
system, but without the hard-limiter in the interrogator. The 
results show that both implementations perform adequately over all 
data rates. Removing the hard-limited improves the tracking 
performance slightly, however at the expense of reduced power 
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Fig. 1.4. Detached Payload Return Link Block Diagram. 











control and hence more crosstalk into the low-rate channel and 
is therefore not recommended. 


The results of our analysis also allow to predict the effect 
of losing the detached payload signal. With the proper choice for 
the Costas loop arm filter bandwidths the loop can still track 
the subcarrier properly. However, if the arm filter bandwidth is less 
than approximately one fourth of the payload interrogator bandwidth 
the lock-points of the loop S-curve are shifted by 90 degrees and 
the channels are interchanged at the detector output. 
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2.0 DESCRIPTION OF UnCsim 

2.1 Introduction 

LinCsim is a software package designed to predict the effect 
of certain signal distortions in the transmitter (those specified 
in the TDRSS User's Guide, [1] Table 3-14 as user constraints) on 
the hit error rate performance of the Shuttle/TDRSS links. 

Analytical simulation has been used as the basis for developing 
a computer simulation model which defines the service performance 
capability of the TDRSS as a functionof the specific subsystem 
characteristics. The model includes the dynamic effects of the TORS/ 
ground terminal links and TDRS/Shutt1e links. The analytical char- 
acterization includes all factors which affect TDRSS services for 
representative operational conditions. This includes (but is not 
limited to) the following: channel linear and nonlinear distortions 

(AM/AM and AM/PM effects, TWT backoff, amplifier limiting), oscillator 
phase noise sources, additive thermal noise sources, filtering effects 
in the TDRS and ground station, demodulation and despreading losses 
for balanced and unbalanced QPSK signals, bit synchronization effects 
and convolutional coding performance. At present the program has the 
capability of predicting the bit error rate performance for S-band 
single access (SSA) and Ku-band single access (KSA) forward and return 
links. Both coherent and noncoherent turnaround implementations, oper- 
ating via the hybrid and cross-support tracking modes are available. 

The analytical simulation program can be used for several purposes: 

(1) To perform a detailed evaluation of several key com- 
munication features of the Shuttle services in order to 
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ensure the consistency of TDRSS and Shuttle user space- 
craft transponder performance specifications. 

(2) To verify that the TDRSS return link and tracking services 
are provided without degradation for a Shuttle user with 
transmitted signal characteristics within the Shuttle 
user constraints of NASA Specifications [1]. 

(3) To simulate the full range of permissible Shuttle space- 
craft characteristics in order to meet the TDRSS achievable 
data rate, 

(4) As a verification tool of Shuttle spacecraft characteristics 
with TDRSS operational. 

In order to fulfill these purposes the program must accurately 
model the true signal format, filtering, nonlinear channel effects 
and any other parameter affecting data transmission or tracking 
performance. Input and output must be in terms of parameters relevant 
to the design engineer. In the following we will discuss the models 
used for the signals and the transmission path and the general set-up 
of the program. Typical results obtained from LinCsim will be given 
in Chapter 3. 

2.2 Link Models 
2.2.1 Introduction 

This chapter presents the models on which the LinCsim computer 
simulation is based. There is a basic channel model which is general 
enough to emulate all forward and return services. This generic link 
model is presented first, followed by channel descriptions for the 
individual services. To complete the channel definition, a list of 
the design values of all parameters associated with the link or one 
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of Its subsystems is included. These design values are based on 
actual system parameter values whenever they are available; other- 
wise, an estimate was used. The source of each value is indicated. 

2.2.2 Channel Model Description 

LinCsim is based on a generic channel model which can be 
modified to fit any particular service (see Sections 2.2.3 and 
2.2.4) by specifying the parameters for each subsystem or by 
bypassing it if necessary. This model is depifcted in Figs. 2.1 to 
2.3. In Fig. 2.1 d.^ and d^ are two binary independent data streams 
(except for the K-band return link where d^ is a subcarrier with 
UQPSK modulation). They may be routed through convolutional encoders 
and NRZ-to-Bi phase converters. Then they may be modulated by the 
independent pseudo-noise sequences PN-j , PN 2 before they modulate the 
IF carrier in the modulator. The modulator output is spectrally 
shaped by the modulator filter. Any linear signal distortions occurring 
in this part of the transmitter are modeled by a distortion filter. 

The signal then passes through a high-power amplifier which distorts 
the signal in a nonlinear fashion. On the uplink, white Gaussian 
noise is added to the signal. 

At the input to the TORS repeater the signal and noise are 
filtered. The power level is then adjusted to the specified value. 

The signal may be combined with other signals in a frequency division 
multiplex format. A nonlinear TWT amplifier provides the necessary 
power for the downlink and an output filter limits the transmitted 
spectrum. At the input to the receiver another white Gaussian noise 
is added. PN-spread channels are then despread before the carrier 
phase and bit timing are extracted. The soft-detected bits are 
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Figure 2,1, Transmitter Model, 
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hard! imi ted for uncoded data or Vi terbi -decoded for coded links. 

2.2.3 Forward Link Models and Design Values 

The models representing the forward service S-band and K-band 
links are shown in Figs. 2.4 and 2.5, respectively, together with 
some of the link characteristics. The design values for the 
parameters defining each one of these models are summarized in 
Tables 2.1 and 2.2. These values are stored in the program and are 
used if the operator does not choose to override them. 

2.2.4 Return Link Models and Design Values 

The models for the return link services and their important 

/ 

characteristics are shown in Figs. 2.6 through 2.8, Two different 
models are defined, one based on a dedicated TWT amplifier (i^e., 
not shared with other signals), the other assuming a shared TWT 
amplifier (composite link model). For the K-band link these models 
correspond to the dedicated and composite links, respectively. For 
the S-band link the dedicated link model is included to provide a 
means of assessing the degradation due to the intermodulation effects 
in the nonlinear amplifier. The two models are discussed in more detail 
in Section 2,2.5, 

The design values for all parameters of the return links are 
summarized in Tables 2.3 and 2.4. Just as for the forward links, 
these are the default parameter values stored in the program. 

2.2.5 Dedicated vs Composite Channel Model for- the S-band Return 
Link 

The S-band dedicated link model is provided as a tool to estimate 
the effects of crosstalk generated in the shared nonlinear TORS 
amplifier on the bit error rate. To allow such a comparison the two 
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Figi. 2,4., S-Band Shuttle Forward Link 
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Table 2. 1 PARAMETER OESTRN VALUES FOR RFR PROGRAM - SHUTTLE FORWARD LINK SSA 


PARAMETER 


DATA RATE 


POWER SPLIT 


MODULATOR PHASE 
IMBALANCE (BPSK) | 

REL. PHASE BE- 
TWEEN I&a CHANNELS 

DATA ASYMMETRY 

PN ASYMMETRY 

DATA SKEW 

PN SKEW 

MODULATOR GAIN 
IMBALANCE 

XTR GAIN FLATNESS 

XTR GAIN SLOPE 

XTR PHASE NON- 
LINEARITY 

XT FILTER BW 

XT FILTER ORDER 

XT FILTER RIPPLE 

XTR AM/ AM 

XTR AM/ PM- •. 

TORS TWT MAX AM^M " 

TDRS TWT AM/ AM 

TORS TWT AM/ PM 

TDRS FRONT END 
FILTER BW 

OKTO 

Rimr 

TDRS TRANSMIT 
•FILTER BW 

mm 

ITTPPLL 

PN TIMING OFFSET 
PN tImING JITTER 
DATA TIMING offset" 


-DATA TIMINGOITTEW 


present 

VALUE 


72 Kb/sec 


+3^ 


WHERE 

OBTAINED 


DATE 

OBTAINED 


OLD VALUE SOURCE OLD VALUE 


Harris PDR 3/30/78 


Harris PDR 3/30/78 


Jz Kb/ sec S-805-1 


1 : 0 


S-805-1 


+ 3' 


S-805-1 


-<•3° peak 


Harris PDR 3/30/78 


+3° peak S-805-1 


+3% 


Harris PDR 3/30/78 


+3% peak S-805-1 


< 1 % 


Harris PDR 3/30/78 


+1% peak S-805-1 


N/A 


< 1% 


Harris PDR 3/30/78 


N/A 


S-805-1 


+0.25 dB 


Harris PDR 3/30/78 


+0.25 dB peak S-805-1 


+0.8 dB peak S-805-1 


.3dB^ MHz 


Harris PDR 3/30/78 


+.ldB/MHz peak S-805-1 


Harris PDR 3/30/78 


+8.6°peak/7MHz S-805-1 


4.5 MHz 


Lin Com estimate 


Lin Com estimate 


0.1 dB 


LinCom estimate 


1 dB / dB 


LinCom estimate 


< 4° dB 


Harris PDR 3/30/78 


LinCom estimate 


10°/dB 


TPM 


3/78 


lOVdB 


LinCom estimate 


0.2 dB/dB I LinCom estimate 


6VdB 


TPM 


3/78 


6°/dB 


LinCom estimate 


20 MHz 


Harris PDR 3/30/78 


20 MHz 


Wu Proposal 7/76 


Lincom estimate 


JLX 


JnCom estimate 


4.5 MHz 


Lin Com estimate 


T 


LinCom estimate 


JLZ. 


I inCnm ps+ima+e 


J 


inCnm 


1 % peak 


S-805-1 


2 % 

1% peak 


LinCom estimate 
T-^05-1 


-IQ- 





























































































2.1. PARAMETER DESIGN VALUES FOR BER PROGRAM - SHUTTLE FORWARD LI NK SSA 


PARAMETER 

PRESENT 

VALUE 

WHERE 

OBTAINED 

DATE 

OBTAINED 

OLD VALUE 

SOURCE 
OLD VALUE 






XTR EIRP: UNCODED 
DATA 

— 

-■- 

IHl 

— 

1. 

XTR EIRP.: COOED 
DATA 

71 . 1 dBW 

Harris PDR 

3/30/78 

7T.1 dB 

WU Proposal , 7 

SPACE LOSS 




208.1 

WU Proposal , 1 

POLARIZATION LOSS 




0.1 dB 

WU Proposal , / 

ATMOSPHERIC LOSS 




0.8 dB 

WU Proposal , 7 

TORS G/T 


• - ' 


10. dB/°K 

WU Proposal , 7 

TORS TWT MAX 
3UTPUT POWER 

14.11 dBW 

TPM 

3/78 

14.1 dBW 

WU Proposal , 7 

FDRS TWT OUTPUT 
BACKOFF 

■ 



0 ■ ■ 

LinCom estimat 

^OWER ALLOCATION j 

0 dBW 

TPM 

3/78 

0 dBW 

LinCom estimat 

:HANNEL POWER 

14.11 dBW 

TPM 

3/78 

14.1 dBW 

LinCom estimat 

FDRS HARDWARE LOSS 



' • 

3.0 dB 

LinCom estimat 

FDRS ANTENNA GAIN 

s. 




35.4 dB 

WU Proposal , 7 / 

’OINTING LOSS 




0.5 

TDRSS Users Go 

iPACE LOSS 

191.4 

EBHIHB 

1/77 

191.6 

TDRSS Users Go 

'OLARIZATION LOSS 




0.5 

TDRSS Users Go 

iTMOSPHERIC LOSS 




0 

SS RF Circuit N 

1X ANTENNA G/T 




- -30.1 

SS RF Circuit N 


— ■ II I 1 

PHASE NOISE BUDGET 


T OSCILLATORS. 1-10 

Hz 



1.5° rms 

S-805-1 

10 to 32 Hz 




— IBI IMI 

S-805-1 

32 Hz to 1 KHz 




ansBiiBn 

S-BUb-i - - 

1 KHz to 6 MHz 





S-805-1 

TORS OSCILLATORS 







X Osnillators 




1.4 

S-805-1 

X CARRIER TRACK- 
NG LOOP - STATIC 

HASP ERROR 

■ 


- 

3° 

LinCom estimate 

BANDWIDTH 




32 Hz 

LinCom estimate 

DAMPING 




2 

LinCom estimate 

SQUARING LOSS 




■ 2 dB 

LinCom estimate 


- 20 - 




































































Figi. 2.5. Shuttle Ku-Band Forward Link 
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Table 2.2. .PARAMETER DESIGN VALUES FOR BER PROGRAM -* .SHUTTLE FORWARD LINK KSA 



PRESENT 

WHERE 

PARAMETER 

VALUE 

OBTAINED 



FA RATE 


WER SPLIT 


DULATOR PHASE 
BALANCE (BPSK) 


L. PHASE BE- 
EEN I&a CHANNEL! 


TA ASYMMETRY 


ASYMMETRY 


TA SKEW 


SKEW 


DULATOR GAIN 
BALANCE 


R GAIN FLATNESS 


R GAIN SLOPE 


R PHASE NON- 
NEARITY 


FILTER BW 


FILTER ORDER 


FILTER RIPPLE 


R M/PM 


RS TWT MAX AMf>M 


RS TWT AM/ AM 


RS TWT AM/ PM ■ 


RS FRONT END 
-TER BW 



RS TRANSMIT 
-TER BW 



TIMING OFFSET 


TIMING JITTER 




PA TIHINGJITTE 




OLD VALUE SOURCE OLD VALUE 


21-6 Kb/sec I SSO ICD/IRN 


1:0 I SSO ICD/IRM 


Harris PDR 


Harris PDR 


Harris PDR 


Harris PDR 


3/30/78 


3/30/78 


3/30/78 


3/30/78 


+3° peak 




+0.75 dB J Harris .PDR 


+:3dB/MHz 


3 


40 MHz 



7°/dB 


7.2°/dB 


0.9 dB 


Harris PDR, 


Harris PDR 


Harris PDR 


N/A 


3/30/78 0.25.dB peak 


0.8 dB peak 


3/30778 O.ldB/MHz 


3/30/78 


3/ 30/ 78 


8 . 6 ° 


50 MHz 



Harris PDR 


TPM 


TPM 


3/30/78 


3/78 


3/78 


r dB/cfB 


37dB 


107dB 


0.2dB/dB 


6°/dB 


50 MHz 


S-805-1 


S-805-1 


S-805-1 


S-805-1 


S.-805-1 


S-805-1 


S-805-1 


S-805-1 


S-805-1 


Lin Com estimate 


Lin Com estimate 


LinCom estimate 


LinCom estimate 


LinCom estimate 


LinCom estimate 


LinCom estimate 


LinCom estimate 


WU Proposal , 7/76 
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2-2 .PARAfCTER DESIGN VALUES FOR BER PROGRAM - SHUTTLE FQRWARr)-| TNI^ kSA 



PARAMETER 

PRESENT 

VALUE 

WHERE • 
OBTAINED 

DATE 

OBTAINED 

OLD VALUE 

SOURCE 
OLD VALUE 


1 LINK BUDGET 



XTR EIRP: UNCODED 
DATA 

75.9 dBW 

Harris PDR 


75.9 dBW 

WU Proposal , 7/ 

XTR EiRP; CODED 
DATA 

— 

— 

— 

— 

— 

SPACE LOSS 




208.4 dB 

WU Proposal , 7/7 

POLARIZATION LOSS 




0.1 dB 

WU Proposal , 7/7 

ATMOSPHERIC LOSS 




0.8 dB 

WU Proposal , 7/7 

TORS G/T 





WU Proposal , 7/7 

TORS TWT MAX 
OUTPUT POWER 

2.1 6 dBW 
Xspec) 

TPM 

3/78 

1.8 dBW 

WU Proposal , 7/7 

TDRS TWT OUTPUT 
BACKOFF 

0 dB 

TPM 

3/78 

0 dBW 

LinCom Estimate 

POWER ALLOCATION 

0 dB 

. TPM 

3/78 

0 dB 

TDRS5 Users Guide 

CHANNEL POWER 

2.16 dBW 

TPM 

3/78 

O.dB 

TDRS -Users Gufde' 

TDRS HARDWARE LOSS 

1.68 dBW 

TPM 

3/78 

1 dB' 

TDRS Users 'Guide 

TDRS ANTENNA GAIN 

51.85dBW- 

TPM 

3/78 

52 . dB 

TDRS Users, Guide 

POINTING LOSS 

■ 



O'. 7. dB 

TDRSS Users Guide 

SPACE LOSS 




207,. 7 dB 

TDRSS Users Guide 

POLARIZATION LOSS 




0.3 dB 

TDRSS Users Guide 

ATMOSPHERIC LOSS 




0 dB 

TDRSS Users Guide 

RX ANTENNA G/T 




■5.3 dB/°K 

TDRSS Users Guide 


PHASE NOISE BUDGET 



XT OSCILLATORS 1 to 

10 Hz 



1.5° rms 

S-805-1 


- 10 to 32 Hz 


. 

wmm\ 

■raBHj 

S-805-1 

oc HZ to j KHz 

Z T T7n 3: ■ 1 

- 

, - - * 



S-8U5-I 




- 





TDRS OSCILLATORS 




3° 

S-805-1 

a 

RX OXCILLATORS 




2.8° 

S-805-1 


RX CARRIER TRACK- 
ING LOOP 




■ 


STATIC PHASE ERROR 
BANDWIDTH • 




3° 

LinCom estimate 




110 Hz 

TRW, 8/77 

. 

DAMPING 




2 

LinCom estimate 

- 

SQUARING LOSS 



2 dB 

LinCom estimate 
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Fig. 2.6., S-Band Shuttle Return Link Model for Data Transmission. 

t 

Convolution Coding: Data Coded in Both Mode 1 & Mode 2 

Data Decoding: Rate 1/3, Constraint Lenoth 7 


Data Rates: Mode 1: 96 Kbps 

Mode 2: 192 Kbps 


PN Code: Not Applicable 


Carrier Modulation: Bi-([) L Symbols PSK Modulate the 

Carrier +^tt/ 2 Radians 0Ji34 

















Table 2.3. PARAMETER DESIGN VALUES FOR BER PROGRAM - SHUTTLE RETURN LINK SSA. 


PARAMETER 

PRESENT 

VALUE 

WHERE 

OBTAINED 

DATE 

OBTAINED 

OLD VALUE 

SOURCE OLD VALUl 

DATA RAtE - I 




192 Kbps 

S-805-1 

DATA RATE - Q 


- - 


0 .. 

S-805-1 

POWER SPLIT 




1:0 

S-805-1 

MODULATOR PHASE 
IMBALANCE (BPSK) 

11° 

JSC 

11/78 

+3° 

S-805-1 

REL. PHASE BE- 
TWEEN I &a. CHANNELS 




N/A 

S-805-1 

DATA ASYMMETRY 

3.8% 

JSC 

WBBSM 

+3% 

S-805-1 

PN ASYMMETRY 




N/A 

S-805-1 

DATA SKEW 




N/A 

S-805-1 

PN SKEW 




N/A 

S-805-1 

MODULATOR GAIN 
IMBALANCE 

.1 dB 

JSC 

n/78 

+0.25 dB 

S-805-1 

XTR GAIN FLATNESS 

.4 dB 

JSC 


+0.3 dB 

S-805-1 

XTR GAIN SLOPE 





S-805-1 

XTR PHASE NON- 
LINEARITY 




+3° 

-w _ 

S-805-1 

XT FILTER BW 

100 MHz 

JSC 

n/78 

1 .2 MHz 

S-805-1 

XT FILTER ORDER 




4 

LinCom estimate 

XT FILTER RIPPLE 




0.1 dB 

LinCom estimate 

XTR AM/AM -■ 




T dB/dB 

LinCom estimate 

XTR AM/ PM 

14°/dB 

JSC 

11/78 

12°/dB 

S-805-r 

TORS TWT MAX AI^M 

r 



10°/dB 

LinCom estimate 

TDRS TWT AM/ AM 


. 


0 dB/dB 

LinCom estimate 

TORS TWT AM/ PM 




6°/dB 

LinCom estimate 

TDRS FRONT END 
FILTER BW 




10 MHz 

LinCom estimate 


■■■■■ 

■■■■Hi 


4 

LinCom estimate 





on dB 


TDRS TRANSMIT 
FILTER BW 




621 MHz 

LinCom estimate 

ORDER 




2 

LinCom estimate 

RTPPTE 




0.2 dB 


PN TIMING OFFSET 




1% 

LinCom estimate 

PN TIMING JIHER 




1% peak 

S-805-1 

..DATA TIMING OFFSET 




1% 

LinCom estimate 

; DATA TIHIN6JITTEF 

6% 

JSC 

n/78 

1% peak 

S-805-1 
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jauic c.j. rn;%MiMLit:K Ut.:>IbP< VALUbb i-UK BtK FKU(iKAN - :aHUilLt KhlUKN LINK 5SA 


PARAMETER 


XTR EI.RF: UNCODED • 
DATA 


XTR EI'Rft CODED 
DATA 


SPACE LOSS 


POLARIZATION LOSS 


ATMOSPHERIC LOSS 


TDRS 6/T 


TDRS TWT MAX 
OUTPUT POWER 


IDRS TWT OUTPUT 
BACKOFF 


POWER ALLOCATION 


CHANNEL POWER 


TDRS HARDWARE LOSS 


TDRS ANTENNA GAIN 


POINTING LOSS 


SPACE LOSS 


POLARIZATION LOSS 


ASPHERIC LOSS 


RX ANTENNA G/T 


PRESENT WHERE DATE SOURCE 

VALUE OBTAINED OBTAINED OLD VALUE OLD VALUE 



OSC, COHERENT 
RNAROUND 


3 Hz to 1 kHz 


kHz to 6 MHz 


OSC NONCOHERENT 
RNAROUND 


to 10 Hz 


3 to 100 Hz 


DO Hz to 1 kHz 


kHz to 6 MHz 


RS OSCILLATORS 


OSCILLATORS 


CARRIER TRACKING 
3P 


TATIC PHASE ERROR 


\NDWIDTH 


WING 


lUARING LOSS 



40.3 dB/°K 


WU Proposal , 7/76 


WU Proposal ,"^7/76 


WU Proposal , 7/76 


WU Proposal , 7/76 


WU Proposal , 7/76 


WU Proposal , 7/76 


WU Proposal , 7/76 


WU Proposal, 7/76 


WU Proposal , 7/76 


PHASE NOISE BUDGET 




DGl; 2 dB 
DG2: 8 dB 


Lin Com Estimate 


TRW, 8/77 


LinCom Estimate 


LinCom Estimate 
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ransmitter. 
e 2 

dth {Predetection) 50 MHz 0^35 



















able 2.4- PARAMETER DF.^^TRN VALUES FOR rfb prorpam . 


'V. (Z. 

SHUTTLE RETURN LINK k~.<^A 




PARAMETER 

PRESENT 

VALUE 

WHERE 

OBTAINED 

DATE 

OBTAINED 

OLD VALUE 

SOURCE OLD VALUE 

\TA RATE - I Ch.‘ 
- Ch.J 
Ch..- 


- - 

- 

192 Kb/s 
2-Mb/s. 

50 Mb/s ■ 

■ S-805-T- 

WER SPLIT 




2a:-4: 1 

S-805-I 

)DULAT0R PHASE 
BALANCE (BPSK) 




+3^* 

S-805-1 

•L. PHASE BE- 
lEEN I &CL CHANNEL 

S 



+3° 

S-805-1 

JA ASYMMETRY 




+3% 

S-805-1 

! ASYMMETRY 




— mr — — 

T-805-1 

TA SKEW 




n/a' 

S-805-1 

SKEW 




N/A . — 

S-805-1 

DULATOR GAIN 
BALANCE 




+0.25 dB 

S-805-1 

R GAIN FLATNESS 




+0.3 dB 

S-805-1 

R GAIN SLOPE 




+0.1dB/MHz 

S-805-1 

R PHASE NON- 
NEARITY 




+3° 

S-805-1 

FILTER BW 



1 


• LinCotn Estimate' 

FILTER ORDER 





LinCom estimate 

FILTER RIPPLE 




O.I 

LinCom estimate 



f 


1 dB/dB 

LinCom estimate 





12°/dB 

S-805-1 

cs IWT MAX Af^M 



- 

lOVdB 

LinCom estimate 

5S TWT AM/ AM 




0 dB/dB 

LinCom estimate 


- . 

■* 


6VdB 

LinCom estimate 





225 MHz 

" ' }\ ■ ■ - 

WU Proposal, 7/76 

BSittiliuSl 



HHII 


Linoom II mite 

IKANSMIT 
TER BW 





LinCom estimate 

UKUfcK 

RTITO 




■SBBSi 


TIMING OFFSET 



1 

i^HI 

i^^Sl 

LinCom estimate 
LinCom estimate 

TIMING JITTER 




1 ^peak 

S-805-1 


* 




LinCom estimate 

A tlMING JITTER 




1% peak 

S-805-1 
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PARAMETER 


PRESENT 

VALUE 


WHERE 

OBTAINED 


DATE 

OBTAINED 


OLD VALUE 


SOURCE 
OLD VALUE 



SPACE LOSS 


POLARIZATION LOSS 


ATMOSPHERIC LOSS 


TORS G/T 


TORS TWT MAX 
OUTPUT POWER 


TORS TWT OUTPUT 
BACKOFF 


POWER ALLOCATION 


CHANNEL POWER 


TORS HARDWARE LOSS 


TORS ANTENNA GAIN 


POINTING LOSS 


SPACE LOSS 


POLARIZATION LOSS 


ATMOSPHERIC LOSS 


RX ANTENNA G/T 


TX OSC, COHERENT 
RNAROUND 


TX OSC, NONCOHERENT 
TURNAROUND 


1 to 10 Hz 


10 to 100 Hz 


100 Hz to 1 kHz 


1 kHz to 150 MHz 


TORS OSCILLATORS 


RX OSCILLATORS 


RX CARRIER TRACKING 
LOOP 


STATIC PHASE ERROR 


BANDWIDTH 


DAMPING 


SQUARING LOSS 



-2.5 dB 


TPM#8 


0 dB 

SSRF CircuitMarc 

22.6 d8/°K 

SSRF Circuit Mar 

12.5 dB 

WU Proposal , 7/ 

Composite, 2. dB 
Dedicated,.5dB 

WU Proposal , 7/ 

-2.4 dB 

WU Proposal , 7/ 



0.7 dB 


207.7 dB 


0.1 dB 


0.8 dB 


40.3 dB 


WU Proposal , 7/ 


WU Proposal , 7/ 


TDRSS Users GuT 


TDRSS Users Gui. 


TDRSS Users Gu1( 


jTDRSS Users Guic 


PHASE NOISE BUDGET 




3 


1000 Hz 


2 


DGl : 2 dB 
DG2; 8 dB 


LinCom Estimati 


TRW, 8/77 


LinCom Estimati 


LinCom Estimat( 
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models are made to agree in every respect except for the sharing 
of the amplifier. In particular, the S-band channel power into the 
TWTA and the amplifier characteristic and maximum TWT output power 
are identical. 

2.2.6 Ideal Channel Model 

Along with the error rate for the specified channel the program 
prints the error probability for an "ideal channel." This is assumed 
to be an infinite bandwidth, linear channel with additive white 
Gaussian uplink and downlink noise as the only distortions. Uplink 
and downlink carrier-to-noise ratios as well as the repeater input 
noise bandwidth {affecting the power robbing in the satellite) are 
matched to the actual channel. The printed error rate is therefore 
the optimum performance achievable with the given link power budget 
and a nonprocessing satellite. 

2.3 Subsystem Models 

2.3.1 Introduction 

This chapter extends the channel model description to the sub- 
system level. The implementations chosen are based on the actual 
Shuttle, TORS or ground station hardware implementation whenever 
such information was available. Otherwise, a reasonable state-of-the- 
art guess was used. The Shuttle transmitter implementation has enough 
flexibility to model the user constraint distortions. 

2.3.2 Modulator 

The modulator accepts two binary switching waveforms and phase- 
modulates them onto a carrier. The phase transitions are instantaneous. 
The resulting waveform may exhibit the following types of distortions: 
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•Gain Imbalance 

•Phase Imbalance {for BPSK) 

•Monorthogonal Channels (for QPSK and UQPSK) 

•Data Asymmetry 
•PN Asymmetry 

•I/Q Data Skew (for equal rates and no PN spreading) 

•I/Q PN Skew 
•Data Bit Jitter 
•PN Chip Jitter 
•Carrier Phase Jitter 
which are defined in Section 2.5.3. 

2.3.3 Filters ' 

All filters are modeled as Chebyshev filters. They are defined 
by the three parameters bandwidth, ripple and number of poles. The 
number of poles refers to the transfer function of the equivalent 
baseband filter. The ripple is defined as 20 log(v.|/v 2 ) where v^V 2 
are shown in Fig. 2.9. The bandwidth is the two-sided ripple-bandwidth 
as shown in the figure, i.e. it does usually not agree with either 
the noise bandwidth or the 3 dB bandwidth. However, if the ripple 
is specified as 0 dB a Butterworth filter characteristic results and 
the 3 dB bandwidth is used. 

In the transmitter some linear distortion may be specified in 
terms of the following parameters: 

• Gain Flatness 

• Gain Slope 
•Phase Nonlinearity 
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The transfer characteristic used is described in Section 2.5.3. 
In order to obtain the phase nonlinearity specified and no more, 
the transmit Chebyshev filter has a linear phase characteristtc. 

All other filters exhibit the unequalized Chebyshev phase char- 
acteristic. 

Note that the filter in the transmitter inherently sets the 
following parameters which are specified as user constraints in 
Ref. [1]: 


•Data Transition Time 

•Minimum 3 dB Bandwidth Prior to Power Amplifier 
2.3.4 Power Amplifiers 
2. 3. 4.1 Introduction 

The power amplifiers are modeled as memory less devices with a 
gain characteristic which may depend on the instantaneous value of the 
signal envelope and with an envelope-dependent phase shift between the 
input and output. The transmitter HPA characteristic therefore sets 
the follov/ing distortion parameters: 

•AM-PM 


•Data Transition Induced PM 

There are three different characteristics built into the 
program: a linear amplifier, a measured TWT characteristic and a 

characteristic with constant AM-AM and AM-PM distortion. These 
models are discussed below. An important parameter for a nonlinear 
power amplifier is the operating point. Its definition is discussed 
at the end of this section. 

2. 3. 4. 2 Linear Amplifier 

The linear amplifier has a constant arbitrary gain and no phase 
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shift. The input power is automatically adjusted so that the output 
power agrees with the specified amplifier output power, viz., the 
output saturation power (dBW) minus the output backoff (dB)‘. It is 
understood that the terms “output saturation power" and "output 
backoff" are meaningless for a linear amplifier; however, they are 
used in the above definition of output power in order to reach a 
closer agreement with the TWT amplifier operating point definition 
which is given below. 

2. 3,4. 3 TWT Characteristic 

The measured TWT characteristic corresponds to the Hughes 261H 
tube and is shown in Fig. 2.10. In order to provide more flexibility 
the model includes two parameters which allow one to adjust the scale 
on the gain and phase nonlinearity separately to meet the output 
saturation power and maximum AM-PM specifications. 

The- AM-PM distortion is defined as 


AM-PM 


dg(R) 

^20 log^Q RJ 


In(lO) Rdg(R) r^egi 
20 dR *■ dB-* 


The normalized AM-PM function is plotted in Fig. 2.11. Note that it 
peaks at approximately 8 dB input backoff. This means that the worst 
degradation due to AM-PM can be expected if the TWT input signal power 
is 8 dB below the saturation power. 
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Figure 2,10. Normalized TWT AM-AM and AM-PM Characteristics. 
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2. 3. 4. 4 Constant m-m, AM-PM Model 
■The AM-AM distortion is defined as 


d[20 log,Qf(R)] 

AM -AM = 

d[20 log^Q R] 

and the AM-PM distortion is 


R df(R) pdB^ 
fW dR ‘-dB-* 


AM-PM = 


dg(R) 


Using 


d[20 log^Q R] 


_ in(lO) Rdg(R) r^ieg-, 

20 HD 1- /<D-1 


dR 


dB- 


f(R) = R' 


AM-AM 


g(R) = 


AM-PM 20 log^Q R R > e 

AM-PM 20 log^Q e R < e 


where e is some small {with respect to the r.m.s. signal) number we 
obtain a characteristic with a constant AM-AM and AM-PM value 
over all operating points. This model is implemented with AM-AM 
and AM-PM as parameters. A typical characteristic is plotted in 
Fig. 2.12. This model can also simulate a hard-limiter by setting 


AM-AM = AM-PM = 0. 


As in the case of the linear amplifier the terms "output 
saturation power" and "output backoff" are not applicable to this 
characteristic. They are however used to define the output power in 
the same way as for the linear amplifier. 
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Figure 2.11. Normalized AM-PM Characteristic. 
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Figure 2.12. Normalized Characteristic for Constant AM-AM, AM-PM. 


•78 0025 




om 


2. 3. 4. 5 Definition of Operating Point 

For the linear amplifier and the constant AM-AM, AM-PM model 
performance is insensitive to the amplifier operating point, viz,, 
the power level of the signal into the amplifier. For the TWT 
model, however, it is an important parameter. Here the operating 
point is characterized by the input backoff, which is defined as 
the ratio of the input power (signal plus noise) to the input 
saturation power, where the input saturation power is the power of a 
CW input signal which produces the maximum output power. The input 
backoff is computed from the specified output backoff on the basis of 
a noise-free CW signal. 

2.3.5 PN Despreader 

The despreader implementation is shown in Fig. 2.13. For a 
four-phase signal, only one quadrature component (the higher-power 
channel for unbalanced QPSK signals) is used to find the PN epoch. 
Both PN sequences are then derived from the same PN generator, 

2.3.6 Carrier Recovery 

The carri er tracking loop is modeled as a Costas loop for BPSK 
signals and for PN-spread signals and as a two-channel Costas loop > 
with hard-limiters (Fig. 2.14) for balanced and unbalanced QPSK 
signals. For PN-spread QPSK or UQPSK signals a single loop tracks 
the carrier in one of the quadrature channels and its phase estimate 
is used for the demodulation of both channels. 

2.3.7 Data Clock Recovery 

The data tracking loop implementation is shown in Fig. 2,15, 

The baseband input signal is hard-limited and multiplied with a 
delayed copy of itself. This produces a line in the spectrum which 
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is being tracked by a phase-locked loop. The delay t can be 
selected; a typical value is half a bit time for NRZ symbols and 
one quarter of a bit time for biphase symbols. The resulting spectral 
line is at the data rate or twice the data rate, respectively. The 
loop bandwidth is set at 1/1000 of the data rate or 10 Hz, whichever 
is more. 

2.3.8 Data Detector 

The data detector is modeled as a hard-limiter for uncoded data 
and as a three-bit quantizer for convolutionally encoded data. 

2.4 Definition and Modelling of User Constraints 

2.4.1 Introduction 

This section lists all the user constraint parameters specified 
in [1], Table 3-14, along with their official definition ([1], 

Appendix I) and the model used in LinCsim. An effort was made to 
find for every parameter the worst-case model meeting the 
specifications. 

2.4.2 User Constraints 
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Definition: The steady state phase difference between the phase 

separation of the BPSK modulated +1 and -1 data bits 
and 180 degrees. 

Model : 



Y = Phase Imbalance 




=£-„(2 


om 


Modulator Gain Imbalance (dB) 

1 . BPSK 

Definition : 10 log P-|/Pg = 20 log[ ^^^^] where P-j is the power in an 

all 1 's data pattern as measured at the modulator output, 
Pq is the power in an all O's data pattern as measured 
at the modulator output, X is the steady state amplitude 
gain of the modulator phase state one, X+aX Is the steady 
state amplitude gain of modulator phase state two, and 
AX > 0. 


Model : 



|20 log(v^/V2)| = modulator gain imbalance 

1 2 
^[max(v, ,v )] = 


design modulator output power 


2. QPSK and UQPSK 
Definition : 20 log 

states have steady state amplitude gains X, X+aX^ , 
X+AX 2 and X+aX^, and where AX^ ^0. 

Model : 


‘X+(AX.)f 

X 


where the four modulator phase 



20 loq(v /v . ) = modulator gain imbalance 
‘ max' min" 

/2 = design modulator output power 
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Data Asymmetry 
Definition; 


1. For an NRZ format signal: 

r iength of long bit - length of short bit- . 

’-length of long bit + length of short bit-* ° 


2. For a bi-0 format signal, data asymmetry applies to both 

the entire bit as well as to each half symbol pulse. 

For data bits, data asymmetry is: 

r iength of long bit - length of short bit- i 
“-length of long bit + length of short bit-* 

For half symbol pulse, data asymmetry is; 

[length of long half symbol pulse - length of 
short half symbol pulse]/[length of long half symbol 
pulse + length of short half symbol pulse] x 100% 


Model : 


Each positive-going pulse transition is advanced by 
aTg/2, each negative going transition is retarded by 
oT^/2 with respect to a perfect clock, where = symbol 
time for NRZ, = 1/2 symbol time for Bi-^ data. Then 

Length of long pulse; 12 ( 1 + 0 ) 

Length of short pulse: T^Cl-a) 

T (l+o)-T (1-a) 

Data Asymmetry = . 1 00 = -1 00 
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Phase Nonlinearity 

Definition : The peak deviation of the phase from the best linear fit 
to the actual phase versus frequency relationship over 
the bandwidth of interest, illustrated as follows: 


‘PHASE 


ACTUAL 

PHASE, 


LINEAR FIT 


A0 (MAXIMUM) 



BANDWIDTH OF lNTEREST-*j ^PREOUENCY 


PHASE NONLINEARITY (£l0) DEFINITION 


The phase nonlinearity is modeled as a sinusoidal ripple 
within the specification bandwidth B. Outside, it flares out 
linearly as shown below. The ripple period is equal to the 
data rate, the amplitude is equal to the specified phase 
nonlinearity. ^ 
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Gain Flatness 

Definition : The peak deviation of the gain from the best horizontal 

fit to the actual gain versus frequency relationship over 
the bandwidth of interest, illustrated as follows: 



Model : If the specified gain slope is sufficiently large the 

logarithmic gain characteristic is modeled as a trapezoidal 
waveform with a peak-to-peak amplitude of twice the 
specified gain flatness AG with a slope equal to the gain 
slope G' specified. Outside of B it flares out. The 
ripple period is equal to the data rate. The resulting 
characteristic is illustrated in (a). If the gain slope is 
too small to obtain the above characteristic a triangular 
waveform is used with the same peak-to-peak amplitude and 
with the given gain slope. The period of this waveform 
is determined by these parameters. The phase is adjusted 
to assure that the positive and negative peaks are reached 

i 

within +B. Outside of B the characteristic flares out 
again, see (b) 
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Gam Slope 

Definition : The rate of change in gain versus frequency 'at any- 


point within the bandwidth of interest, illustrated 
&s follows: 



logarithmic gain characteristic is modeled as a trapezoidal 
waveform with a peak-to-peak amplitude of twice the specified 
gain flatness AG with a slope equal to the gain slope G' 
specified. Outside of B it flares out. The ripple period 
is equal to the data rate. The resulting characteristic 
is illustrated in (a).- If the gain slope is too small 
to obtain the above characteristic a triangular waveform 
is used with the same peak-to-peak amplitude and with the 
given gain slope. The period of this waveform is determined 
by these parameters. The phase is adjusted to assure that 
the positive and negative peaks are reached within +B. If 
the gain slope specified is too small, i.e., G'B < aG the 
slope is readjusted to G‘ = AG/B. Outside of B the 


characteristic flares out again, see (b). 




Definition : AM/PM = worst case - over the range of operating 

points where = RF power input to power amplitier, 
in dB and = RF phase output from power amplifier 
in degrees. 



where 

R = RF input envelope (V) 

e = some small voltage (with respect to r.m.s. 
signal ) 

This results in a constant value for the AM/PM. 
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Data Transition Induced PM 

Definition : Phase variation at the power amplifier output which 

is correlated with the modulation on the carrier. This 
phase variation results from the combined effects of 
band! imi ting followed by AM/PM conversion. 

Model : The data transition induce'd PM is determined by the choice 

of type and bandwidth of the pulse shaping filter together 
with the selected value of transmitter AM/PM. 
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PN Chip Jitter 

Definition : Let eT be the error in seconds between the actual 

chip transition time and the correct transition, when 
T is the chip period and e is a random variable. Then 
PN chip jitter is the square root of the variance of e 
multiplied by 100 percent. 

PN chip jitter is measured as the rms phase error in 
a baseband bandwidth of 1 MHz using a high gain 
second-order PN Code Tracking Loop with a damping 
factor of 0.707 and a two-sided noise bandwidth 
(2Bj^) of 0.2 Hz, This user constraint is based upon 
zero TDRSS forward link PN jitter. 


INPUT 
SJGNAL ^ 

PN CODE 
TRACKING LOOP 
2B^»0.2HZ 

PHASE DETECTOR 
OUTPUT 

FILTER 

MEASURE 
PHASE ERROR 



BW=IMHZ 



Model : This parameter describes the random variation of the PN 

clock phase. A Gaussian distribution is assumed and the 
rms value is set to one third of the peak value specified. 


cJCinC^m. 
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Data Bit Jitter 

Definition : User spacecraft peak clock frequency jitter and 

peak jitter rate (sinusoidal or 3o random) as 
percent of the symbol clock rate. 

Model : This parameter describes the random variation of the 

data clock phase. A Gaussian distribution is assumed. 
Note that this definition deviates from the definition 
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Spurious PM 

Definition : Consider a carrier 

sin[oi^t+0(t)] 

0(t) = PM spurs + phase noise. 

-PM spurs are discrete phase modulated 
components on the carrier, 

-Phase noise is random phase modulation. 
Then each PM spur can be written as 
sin(u^t) 

and the resulting rms spurious PM will be 


where the sum is taken over all 1 for which w. 
is within the specified frequency range. 

Spurious PM can be determined by measuring the power 
in each spur and calculating the rms value for these 
powers . 

Model : This parameter is not modelled presently. 


B. 


Spurious PM = 
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Definition : 


Model : 


The short-term EIRP stability is determined by 
measuring the peak-to-peak variations in signal 
EIRP in dB and dividing by the specified averaging 
periods in seconds. 

This parameter has no effect on the bit error performance 
and is therefore not modelled. 




Spurious Outputs 

Definition : The sum of the power in all spurs within the 

specified bandwidth measured relative to the total 
signal power. (dBc indicates dB below total signal 
power). Data bandwidth is the bandwidth between the 
first nulls of the signal. That is, two times the 
symbol rate for NRZ data and four times the symbol 
rate for Bi-d formatted data. The channel bandwidth 
is 6 MHz for MA; 10 MHz for SSA, S-band Shuttle, and 
S-band IF service; and 225 MHz for KSA, Ku-band 
Shuttle, and Ku~band IF service. 

Model : This parameter is not modelled presently. 



Frequency Stability (Peak) 

» 

Definition ; The peak instantaneous carrier frequency deviation 
from the nominal carrier frequency (Af) normalized 
to the nominal carrier frequency as observed over 
the specified time interval of interest. This 
includes frequency deviation due to all sources 
including deviations induced by environmental 
changes. (This parameter only applies, to the non- 
coherent modes of operation.) 

Model: This parameter has no effect on the bit error performance 

•9 

and is therefore not modelled. 
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Incidental AM 

Definition : Consider an ideal RF signal of the form: 

A cos[(o^t+?i(t)] 

Then a signal with amplitude modulation m^- at 
frequencies with phase will be of the form: 

A[1 + m. cos(oi.t+0. )]cos[&3 t+0(t}] 

1 1 c 

The peak of the amplitude modulation will then occur 
when all of the incidental AM components line up. 
Therefore , the peak incidental AM will be: 

2 X 100% 

Model : This parameter has little effect on the bit error performance 

and is therefore not modelled. 
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Axial Ratio 

Definition : For circularly polarized antennas, the electrical field 

vector usually produced describes an ellipse instead of 
a circle. The axial ratio is a measure of ellipticity and 
is the ratio of the minor axis of the ellipse to the major 
axis. 




Phase Noise 

Definition : Phase noise is the random component of th.e 'phase 

modulation on the carrier. If the spectrum of the 
phase noise is taken as S^{f), then the rms phase 
from f^ to f2 is: 

Phase Noise = 



Phase nosie can be determined by measuring the power 
in each frequency band and then subtracting the power 
in the spurs within that band. 


Phase noise values for the coherent turnaround mode 
assume no phase noise on the signal received by the 
user and therefore represent the phase noise added 
by the user. Values indicated for the noncoherent 
mode represent total output phase noise of the user. 
Model : The phase noise is modelled as a Gaussian random process 

with a flat spectrum in -each of the frequency bands 
specified. 
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I/Q Data Skew 

Definition : Relative time delay of the data transitions between 

I and Q channels measured as a percent of a bit time. 
This parameter is determined by data, clock, and 
modulator characteristics. 

Model : This parameter applies to balanced staggered QPSK only. 

It defines the relative time delay between the pulse center 
in the one channel and the pulse transition in the other 
channel. Due to the ground station hardware implementation 
this parameter has a negligible effect on the performance 
of an unspread link and no effect on a PN spread link. 
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Permissible EIRP Variation 

Definition : Range over which the user EIRP measured along the 

user/TDRS 1 ine-of-sight may vary without requiring 
user reconfiguration.' Performance is determined from 
user transmitter power variation, transmit antenna 
pattern, worst case user orientation, and maximum 
variation in range between the user and TORS over the 
duration of a pass. 

Model : This parameter has no effect on the bit error performance 

and is therefore not modelled. 
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Frequency Error of 8.5 MMz Subcarrier 

Def initi on : Difference betv/een subcarrier frequency and 8.5 MHz. 

Model : The frequency error itself has no effect on the bit error 

performance and is therefore not modelled. The effect of 
the loop stress in the receiver is modelled through the 
resulting static phase offset. 
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Subcarrier Phase Noise 

Definition: Phase noise on the 8,5 MHz subcarrier (see phase 


Model : 


noise definition). 

The subcarrier phase noise is modelled as a Gaussian random 
process whose power is concentrated outside the subcarrier 
tracking loop bandwidth but considerably below the data 
rate (worst-case assumption). 
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I/Q PN Skew 


Definition ; Relative time delay of the chip transitions between 

I and Q channels measured as a percent of a chip time. 
Model : Relative time delay between the chip center in one channel 

and the chip transition in the other channel of a SQPN 
signal . 


cJCinC^ 


m. 
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PN Asymmetry 

Definition : PN chip asymmetry is defined as: 

r iength of long chip - length of short chip- , , 
‘■length of long chip + length of short chip-' ^ 


Model : 


Each positive going chip transition is advanced by aT^/2, 
each negative going transition is retarded by «T^/2 with 
respect to a perfect clock, where = chip time. Then 


PN Asymmetry 


t^(l+a)+T^(l-a) 


a-100. 
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3.0 LinCsim PERFORMANCE PREDICTION 

3.1 Introduction 

The Shuttle/TDRSS link simulation software package (LinCsim) 
has been used to predict the Shuttle/TDRSS S-band 
return link bit error rate performance based on current link budgets 
user constraint values and TDRS/ground station hardware data. The 
results of this effort are presented in this section. 

3.2 Performance Prediction for Shuttle S-band Return Link 

This section presents the sensitivity of the Shuttle S-band 
return link BER performance to variations of the user constraint 
values based on current link budgets, nominal user constraint value 
estimates and TDRS/ground station hardware data. This simulation 
is based on the link models presented in Chapter 2 and on the system 
parameters listed in Table 2.3. 

The link characteristics assumed are summarized in Table 3.1. 
The nominal values of the user constraints are listed in Table 3.2. 
They agree with the expected Shuttle values obtained from JSC. The 
link budget is reproduced in Table 3.3. It is based on the budgets 
contained in Ref. 1 with some updates obtained from Dr. Kwei Tu. 
Table 3.4 lists the user constraints whose effect on performance was 
studied. 

The results are shown in two different forms. The BER 
curves show the bit error probability as a function of Shuttle- 
to-TDRS link carrier-to-noise ratio (CNR) variation around the 
nominal link budget of Table 3.3. The horizontal distance between 
the bit error rate curve and the design point shown represents 




'cHinC^c 


om 


Tables.!. Shuttle S-band Return Link Characteristics 
Used for LinCsim 


Data Rate 
Data Format 
Coding 

Carrier 

Design Error Rate 


192 Kbps 

BPSK, Biphase, Unspread 

Rate 1/3 Convolutional 
Code 

Noncoherent with Forward 
Link 
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Table3.2. Nominal User Constraint 
LinCsim 

Data Bit Jitter {3a) 

Modulator Phase Imbalance 
Modulator Gain Imbalance 
Data Asymmetry 
Phase Nonlinearity 
Gain Flatness 
Gain Slope 
AM/ PM 

3 dB Bandwidth 
Phase Noise 

1 Hz - 10 Hz 
10 Hz - 100 Hz 


Values Used for 

. 6 % 

ir 

.1 dB 
3.8% 

3“ 

,4 dB peak 
.1 dB/MHz 
14°/dB 
100 MHz 

0 *^ 

1 ° 


0 
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Table 3.3. Shuttle S-Band Return Link Power Budget, 


Shuttle-to-TDRS Link 


Shuttle EIRP 


Space Loss 
■Polarization Loss 
TDRS G/T, 


16.7 dBW 
192.1 dB 


9.55 dB/°K 


TDRS-to-Ground Link 

TDRS TWT max Output Power 
TUT Output Backoff 
SSA Power Allocation 
TDRS Hardware Losses 
TDRS ' Antenna Gain 
Pointing Loss 
Space Loss 
Polarization Loss 
Atmospheric Loss 
Ground Station G/T 


13.4 dBW 


-10 dB 
2.46 dB 


45.9 dB 


.65 dB 
207.7 dB 


.1 dB 


.8 dB 

40.3 dB/°K 




Table 3.4. parameters Studied 

Modulator Gain Imbalance 

Modulator Phase Imbalance 

Data Asymmetry 

Data Bit Jitter 

Data Static Timing Offset 

XTR m/m 

XTR AM/PM 

Static Phase Error 
XTR Gain Slope 
XTR Gain Flatness 
XTR Phase Nonlinearity 
Phase Noise R.M.S . ' 




the margin in carrier-to-noise ratio which can be allocated to the 
various subsystems for hardware degradations. The sensitivity curves 
show the increase in Shuttle-to-TDRSS CNR needed to offset the 
performance degradation (relative to the nominal performance) due 
to the variation of a single parameter at the design error rate of 


The error rate curve in Fig. 3.1 represents the BER performance 
of a Shuttle transponder transmitting a perfect signal (i.e. all 
user constraints are set to zero). The CNR loss shown (.9 dB) can 
be attributed to the TORS and the ground station. This BER curve 
is reproduced on all other BER plots and labeled "Perfect User." 

The horizontal distance between this curve and one of the other 
BER curves represents the CNR loss due to the combined effect of 
all the user constraints. For the nominal conditions this loss 
amounts to .75 dB. 

Fig. 3.20 shows that by far the biggest contribution of degradation 
comes from the phase noise which is specified as 10° rms in the 0 to 
270 Hz frequency range. As a worst-case assumption this noise power 
was concentrated in the 100-270 Hz range (i.e. outside the tracking 
loop bandwidth). A more accurate performance prediction could be 
made if the noise power inside and outside the tracking loop 
bandwidth were known separately. 
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Figure 3,4.BER Plot for Phase Imbalance. 
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Figure 3..8. BER Plot for Bit Jitter. 
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Figure 3.12. BER Plot for AM/PM 
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Figure 3.16, BER Plot for Gain Flatness. 
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4.0 ATTACHED SPACELAB PAYLOAD BIT JITTER PERFORMANCE 

4.1 Introduction 

In this section the problem of bit jitter on the attached 
spacelab payload link is studied and the effect on the bit error 
rate performance is predicted. The analysis also includes the effect 
of data asymmetry. The results show that the limiting level of 
acceptable clock jitter at the output of the attached spacelab 
payload is determined by the sensitivity of the Viterbi decoder 
to imperfect pulses at the ground station receiver. In Section 4.2 
we introduce the mathematical model we have adopted for our analysis. 
We also present in this section the model we have selected for 
clock jitter, bit jitter and data asymmetry . Based on laboratory 
measurements, we have 'determined the shape of the possible data 
sequences out of the attached spacelab into the Ku-band Space 
Shuttle equipment. These various pulse shapes are also presented 
in Section 4.2. In Section 4.3 we derive, an expression for the 
bit error rate at the spacelab Ku-band transmitter interface due 
to imperfect timing. We then show that for the available SNR the 
effect of an acceptable level of clock jitter at the ground 
receiver has a negligible impact on the overall link performance 
even for a biased sampling time of +^12.5% in the Shuttle repeater. 
(This figure for the maximum bias was given to LinCom by JSC 
personnel. Note, however, that this does not agree with the 
description in [3]. According to this source the worst-case 
sampling point is 5 nsec from the transition. At the highest 
data rate of 50 Mbits/sec this is equivalent to 25% of the bit- 
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time. Hence the maximum bias would also be 25%. ) The coded link 
bit jitter performance is analyzed in Section 4.4. We show in 
this section how a completely soft Viterbi decoder significantly 
ameliorates the effect of imperfect pulses at the output of the 
Ku-band transmitter. Finally, in Section 4.5 we summarize the 
results and present some conclusions on the effect of clock and 
bit jitter on the overall system performance. 

4.2 System Model, Clock Jitter and Data Asymmetry 

The mathematical system model we have adopted for our analysis 
is as shown in Figure 4.1. Notice that although this is a simplified 
model of the overall system it allows us to accurately represent 
the various sources of link performance degradation. It is assumed 
first that the data flow from the attached spacelab payload to the 
Ku-band Shuttle transmitter is at baseband.. The shape of the pulses 
carrying the data has been determined from laboratory measurement (4.2) 
and for the purpose of our analysis they will be assumed to be as 
shown for a typical case in Figure 4,3a Other typical pulse sequences 
are shown in Figure 4.3b. The clock information entering into the 
Ku-band transmitter is assumed to be emanating from an imperfect oscillator. 
The phase detector output e(t) of the clock regenerator circuit will 
then be a random process whose statistics depend on the oscillator- 
phase noise characteristic. If we now sample this process e(t) at 
times t* and t*+T where T is the inverse of the symbol rate, then we 
obtain two random variables e(t*) and s(t*+T). Pictorially we have 
the following situation: 
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Figure 4.1. | MATHEMATICAL SYSTEM BLOCK DIAGRAM I 
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(b) Pulse shape characteristic. Trailing edge. 
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Figure 4.3a. NORMALIZED PULSE SHAPE CHARACTERISTIC 
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We shall assume that e(t*) and e(t*+T) are jointly Gaussian correlated 
random variables and such that 


= Var[e(t*)] = Var[e(t*+T)] (1) 


We now define two new random varibles 


'ck 


= £(t*)-t-e(t*+T)+T 

2T 


and 


^ _ E(t*+T)-£{t*)+I 

Tb = 


( 2 ) 

(3) 


It is clear then that the and are two independent Gaussian 
random variables 


Tck" <'» 


we 


refer in what follows to 



( 6 ) 


as clock jitter and 


inborn 


as bit jitte? and where y is the correlation coefficient between s{t) 
and e{t*+T). Notice that sampling instant of the imperfect pulse 
entering to the Ku-band transmitter from the attached spacelab payload 
will be governed by the random variable 

The output of the convolutional encoder is fed into a driver 
which we assume introduces data asymmetry. The model for data 
asymmetry we have adopted for our analysis is described in.[l] and 
is such that a positive going transition occurs early and a negative 
going transition occurs late relative to the nominal transition time 
Instants. Data asymmetry is then defined as the difference in 
length between the shortest and longest pulses in the sequence 
divided by their sum. For our model we get then 


_ T{1+2(5)-T(1-26) _ 

^ ■ lTi+26 )+tii-25T ■ 


( 8 ) 


where 6 represents the fractional (relative to the nominal symbol 
duration T) increase in positive pulse length due to a single adjacent 
negative pulse- 

4.3 bER at the Spacelab Ku-Band Transmitter Interface 

For this part of the analysis we can concentrate on the mathe- 
matical model of Figure 4.4. The input NRZ pulses into the Ku-band 

transmitter are sampled at random instants of time governed by the 

_2 2 

random variable t^ where t^ = ^ (T ince the clock signal 

is synchronous with the NRZ data, for the purpose of our analysis we 
can assume that, the NRZ pulses are of perfect duration T. The 
probability of error at A conditioned on a sampling instant t^ and 
an input sequence S is 
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P^(t.S) = |erfcj ■ { (9) 

2 

where cr^ is the power of the noise entering into the Ku-band Shuttle 
transmitter from the attached spacelab payload p(t,S) is the pulse shape 
input function which at any given instant of time is a function of at most 
three consecutive data bits. The relationship between pCt^jS) and 
P(t,S) as given in Figure 4.3 is 

P(a,S) = 


We can obtain the unconditional bit error probability by averaging (9) 
over all possible sequences S of length three and the random variable 



where p 


2 


is the signal to noise ratio at the interface and given by 



( 12 ) 


and 

When the sampling instant is biased by some "^bi as ’ then T^j,. has 
to be replaced in (11) by Tck‘'’T5-{as ’ ^bias normalized 
to be between 0 and 1. The relationship in (11) has been programned 
and the result of this computation is shown in Figure '4.5. 



_1 7 I7_ 



BIT ERROR RATE 


Figure 4.. 5. 


BER PERFORMANCE AT KU BAND ENCODER INPUT (A) 





CLOCK REGENERATOR OUTPUT JITTER (%) 



JlinCc 


om 


4"?4 Coded Link Bit Jitter Performance 

We are interested now in the bit jitter performance between 
points C and D of Figured. We shall assume here without any 
significant loss of accuracy in the final answer that the receiver 
is implemented as an ideal integrate and dump detector followed by 
a completely soft Viterbi decoder. We shall also assume that the 
timing at the receiver is perfect. 

Lets assume now that a codeword £= (c^ ,C 2 ,. . • where c^ = jjl 
is transmitted out of the Space Shuttle Ku-band transmitter. We model 
the channel as an ideal AWGN channel, we then have thatthe decoder input 
signal 

r. = x(c) + £ (13) 

where 


r. = 

(r.j,r2»««* 

(14) 

(x-j (c) 

>^2 ^ ’ * • * ^ 

(15) 

£ = 

(n-j ,n 2 s • • • 

(16) 


lets define now for simplicity of notation 

t.j = eC(j-l)T! (17) 

t^j = eDT] + T ( 18 ) 

then it is possible to show after some algebraic manipulation that 
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*u(T-t^-Cj.5T) + Cj.jT[l-u(T-t^j-(j5T)] 


- Cjgj-1 )T+t^.j+Cj_-|6T]u(t^j+Cj_^6T) - Cj(j-1)T 
.D-..(t..+c._^6T)] + c.,T(T-t,.-c.5T)u(T-t,.-Cj6T)} 


If we now define as we did before 




Tp> = 

then from (17), (18) ,(20) & (21) it is possible to argue that t| 
is a slow varying random variable and such that 


jii) = jM = , 

ck 'ck ck 


r(j) -i 


and that ' is a fast varying random variable such that 


E{T^T^} = E{T^}E{T^} for £ j (23) 


We then set 


. = T 


This allows us to rewrite (19) as follows 


cJ^inC^i 
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* "j+l'’-"ck-'bj-'=o«>"<i-^ck-'bj-^/> : 

’‘j*‘^j-l’'^j’VT''ck’''bj^ (23) 

In what follows and 1n order to avoid any boundary condition problems 


we let 


'^O " ‘^N+l “ ^ 


Lets assume now the decoder picks the codeword ^ of the dimension N and c 
is the correct one, then the probability of an error even is 

n / 1 4 . - lli(c)-x(c)l| . 

^d ■ 2 27RJ" 

I 

where d is the Hamming distance between £ and £, 

^ ' ^"^bl • ’■^bN^ 
and |lx.( £)-£(£) II is the Euclidean norm given by 


lkL(£)-x(c)lf = 


> 1 

Zi ’'=J ’=j+l ■■"ck’''bj''''j('=j-l ’’^j+1 ’■'ck’’"bj’^‘ 

( 29 ) 


We now let the correct code word be the all Ts code word, then 


cJ^inC^i 
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so that (24) reduces to 


(31) 

We can now rewrite (27) as 




2 /s 

^ (32). 


In order to obtain the unconditional probability of an error event 


we need to average (32) over and as 




(33) 


The conditional expectation in (33) is in general very hard to find 
whenever N is larger than three. We notice on the other hand that 
the expression inside the square root in (32) is a sum of independent 
random variables and for the K=7 rate 1/2 code we are interested in 
N is always larger or equal to 18. This allows us to use a central 
limit argument to write 




- [y-n ( T|.|, ) f /2a^ ( T_ ^ ) 


ck^ 


/2ira^(T^j^)^ 


f(y)dy (34) 


where 


>-<'ck> = E (35) 


j=l 



_i 9n_ 
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j=l 

j=l . 


f(y) = ^erfc‘,j^/r| 


(36) 


(37) 


The bit error can be found from the probability of an error 
event as given in (33) and the transfer function bound [2], This 
approach is very tedious and usually computationally. impractical . 

We can then get a very good approximation to the bit error probability 
by computing- first the probability of an error event for all those 
paths which are at distance correct path. Lets assume 


there are n such paths, we then find p| such that 


F S '’d. = I erfc{/a7 p*} 

1=1 


(38) 


or in terms of an equivalent bit energy to noise spectral density ratio 


^erfc-' I £ P. 

I" W 'Ifi 


“b ‘ 


W 


(39) 


We then set the bit error probability to be the one that corresponds 
to the decoding of a jitter and data asymmetry free signal transmitted 
over a channel with an £jj/Nq = p^. This approach has been progranmed 
and the result of this computations are shown in Figures 4.6, 4.7, and 4.8. 
The clock and bit jitter parameters in these curves are defined with 
respect to the Shuttle Ku-band transmitter clock regenerator phase 
detector output and the cycle duration time T of the clock 


signal out of the attached space! ab payload. 

*In what follows we let d^ =d... 

tree f 
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Figure 4.8, Sensitivity Curve for Clock Jitter. 
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4.5 Summary and Conclusions 

We have analyzed here the bit and clock jitter performance 
of the attached spacelab payload. We have first shown that the 
effect of an imperfect sampler at the spacelab Ku-band trans- 
mitter interface does not impose any measurable loss in performance 
for rms clock jitter levels less than 10%. On the other hand we 
have also shown that for a jitter-free link the required Ej^/Ng 
tomaintaina BER of 10”^ with 10% of data asymmetry is 4.72 dB. 

This implies a aCNR degradation of 0.52 dB with respect to the 

4.2 dB of Ej^/Ng required for an ideal coded channel to achieve a 
-5 

BER of 10 . In Table 1 we show the additional ACNR (with resepct 

to 4.72 dB) required due to the presence of jitter*. 


Table 1. 


ACNR 

Bit Jitter 

dB 

% 

0.1 

5.70 

0.2 

8.01 

0.3 

9.75 

0.4 

11 .10 

0.5 

12.48 


It is very important to observe that for the levels of the 

jitter shown in Table 1, the BER at the SL-Ku-band transmitter 

-9 

interface, remains below 10 . Its effect on the BER at the 

above mentioned system interface can therefore be ignored. 
Appendix 

In this appendix we attempt to predict the SL HDM return 


*The numbers shown in Table 1 correspond to worst case jitter for 
both eooch and duration .iitter. 
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link performance based on the clock pulse shape characteristic 
and the spectrum characteristic shown in Figure 4.2, which has 
been provided to us by Sid Novosad. From Figure 4.2 we have 
estimated* that the C/Ng present at the HRM clock output is approxi- 
mately 


C/Ng = 104 dB-Hz. 


For this level of C/Ng we have concluded* that the jitter observed 
in Figure 4.2b originates mostly from thermal effects. If we now 
assume that the Shuttle clock regenerator loop bandwidth is 10 KHz, 
it is possible to predict that the jitter at the scrubber output, 
including non-thermal effects, will be less than 2%. For this jitter 
level, the overall performance degradation taking into account the 
0.52 ACNR due to 10% of data asymmetry will be 

aCNR = 0.52 dB + 0.01 dB 

Data Jitter 

Asymmetry 

= 0.53 dB 


*Within the accuracy that this picture allows. 
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5.0 SHUTTLE PAYLOAD BENT-PIPE MODE SUBCARRIER RECOVERY 

5.1 Introduction 

In the bent-pipe mode of Shuttle payload data transmission, 
the data are not detected at the Shuttle payload interrogator. 

They go through the Shuttle repeater, where they are low-pass 
filtered, hard-limited, and modulated onto the high-power channel 
of the subcarrier UQPSK modulatpr. In this section the performance 
of the Costas loop recovering the subcarrier in the ground station 
receiver is studied for this case and compared with an implementation 
which does not include a hard limiter in the bent-pipe link. The 
results show a slight degradation of the loop performance when a 
hard-limiter is used. This advantage is offset, however, by the 
tight power control afforded by the limiter which reduces the 
effects of interference on the other two data channels. 

5.2 Problem Statement 

In the detached-payload bent-pipe mode the binary data stream 
is received by the Shuttle payload interrogator immersed in 
additive white Gaussian noise. This signal goes through the Shuttle 
repeater, where it is low-pass filtered and hard-limited. Then it 
is modulated onto the high-power phase of the 8.5 MHz subcarrier. 

This is shown in Figure 5.1. 


For the analysis of the 8.5 MHz subcarrier recovery, we may 
assume that the carrier has been recovered perfectly. This leads 
to the link model of Figure 5.2. 

We wish to characterize the phase-recovery performance of the 
Costas loop. The theory is developed for an arbitrary signal format 
at the Shuttle repeater output. The results are then applied to 



- 128 - 



-129- 



Figure 5.1 . Three-Channel Transmitter in Bent-Pipe Mode 
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Figure 5.2. Simplified Link Diagram. 
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the case of baseband data with two different implementations for 
the Shuttle repeater: a hard-limiter preceded by a wide low-pass 

filter which does not affect the data signal and an arbitrary low- 
pass filter without hard-limiter. 

5.3 Costas Loop Model (Following Reference 1) 

The input signal to the Costas loop is the random process (see 
Fig. 5.2), 

w{t) = H(p)(s^.(t) + Sq{t) + u(t)) (1) 

where H(p) is the ground station receiver filter, s^(t) is the 
high-power bent-pipe signal, s^{t) is the low-power digital voice 
signal and u(t) is the downlink thermal noise. This signal can 
be written 

w{t) = /Fa(t) sin *(t) + v7b(t) cos ^{t) (2) 

for some random processes a(t) and b(t), where $(t) =(U 2 t+ 0 (t) 
and 0 {t) = J 2 Qt+ 0 Q is the input phase to be estimated. Let $(t) be 

the loop estimate of $(t). In the upper arm of the loop, w(t) is 
multiplied by cos o(t), in the lower arm by sin <l(t). The 
multiplier units have gain /i^ and are insensitive to the double- 
frequency terms. Each output passes through a filter corresponding 
to G(p), yielding the upper-arm signal z (t) and the lower-arm 
signal z^(t). 

z^(t) = v'K.jK^ G(p)[a(t) sin«;>(t) + b(t)cos<f (t)] (3) 

^ /K^ G(p)[a(t)cos f (t) - b(t)sin <^(t)] (4) 
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where ‘/'(t) e $(t)-4>(t) 1s the loop phase error. The two signals are 
multiplied together, with unit gain, to produce the dynamic error 


signal 


z(t) = z^(t)z^(t) 


We assume that‘f(t) varies much more slowly than a(t) and b{t) 
and that the G(p) filter Is sufficiently wide so that 

z(t) - {[(G(p)a(t))2-(B(p)b(t))2] I sin(2<f(t)) 

-EG{p)a(t)G(p)b(t)]cos(2<#) (6) 

The instantaneous frequency of the VCO output is related to z(t) by 

/V 

- K/(p)z(t) + (7) 

where Is a gain constant. Hence, the stochastic integro-differential 
equation of operation of the loop is 

2 - 2KyF{p)z(t) (8) 

Conditioned on <f, z(t) can be partitioned into a nonrandom part 
and a zero-mean random process: 


where 


:(t) = S(<^) + n (t,c^) 


S{c|>) = E{z(t)l^P} 

n,(t,g>) = z(t) - S («f) 


This allows us to rewrite (8) as 

2 = 2fi. - 2K,,F(p)S[g>(t)] - 2KyF(p)n_(t,q>) 
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This equation describes a non-Markovian diffusion process. However, 
under suitable conditions (in particular, if the process n^(t,</) 
is considerably faster than the process <p{t)) it can be approximated 
by a Markov process and Fokker-Planck techniques can be applied to 
characterize the stationary distribution of the modulo-2Tr- reduced 
phase noise process as well as the cycle-slipping rate. For this 
analysis the S-curve S{<P) and the statistics of the equivalent noise 
process n^(t,(^) are needed. 

5.4 Costas Loop S-Curve 

We will assume that JiQ=0 and that the baseband equivalent of the 
H(p) filter is symmetric. We first obtain the input signal to the 
Costas loop. 


We will denote a narrowband signal a(t) with center frequency 

"o by 

^ -i,.. -f- 


d - 



(13) 


a{t) 

= /2 Re[a' (t)e ^ ] 

where a'(t) is 

the baseband equivalent to a(t). 


Defining 

the payload 

signal as 



s^.(t) 

= y(t) sin(oj 2 ^^ 

(14) 

where y(t) is 

normalized 

such that E{y^(t)} = 1 , we find 




s:(t) = -j/T^y(t) 

(15) 



xl(t) H H'(p)si(t) 

(16) 


where H'^{p) is the baseband equivalent of H(p). 
For the low-rpower signal we have 
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Sq(t) = -y^ d^(t)cosU2t) 

(17) 


s'(t) = VT^d^(t) 

(18) 


x'(t) = H‘(p)s^{t) 

(19) 

Now we consider the noise u(t), a narrowband Gaussian process 
with center frequency 1^2 - 


u(t) - i/^u^ {t)cos({ 02 t) - / 2 u 2 (t)sin(u 2 t) 

; (20) 

The 

spectral density, of u(t) is 



S,(f) = 

(21) 

for 

some real function where f 2 = 1 ^ 2 ^ {2i,). We define 



Ry (t) = E[u^.(t)u^.(t+x)] , i = 1,2 

(22) 


R (t) = E[u.(t)u. (t+x)] , i k 
1 k ' ^ 

(23) 

Letting R^, be the inverse Fourier transform of we have 



R, = R„ = Re R , 
U-[ U2 u 

(24) 


= -R , = Im R , 

U 12 U 21 u 

(25) 

We assume' that 1s symmetric about 0, so that 




(26) 


^,2 " ■'’''21 ‘ “ 

(27) 
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Since we assumed that the H’(p) filter is symmetric then H’(p)u-j(t) 
and H'{p)u 2 (t) are real -valued processes. 

Now we go into the Costas loop. 

In the upper arm the input signal is multiplied by i^K^K^cos {o) 2 t-^) 
and the double-frequency terms are dropped. The same result is 
obtained if the baseband equivalent of the signal is multiplied by 

e*^ and the real part is taken. In the lower arm of the loop, 
the input signal is multiplied by /ZK-jK^' sin{o) 2 t-<^) and the double- 
frequency terms are dropped. This is the same as multiplying the 
baseband equivalent by j/K-jK^ e and taking the real part, or 
multiplying it by and taking the imaginary part. 

Since and are each linearly related to s^. , s^, and u, we 
may wri te 


= Zj,(t;s^) + z^(t;s^) + z^(t;u) (28) 

z^(t) = Zs(t;s^.) + z^(t;Sq) + z^(t;u) (29) 

where, for example, z (t;s.) is z (t) when the loop input is just 

w I V* 

H(p)s^.(t}.' 

For s^. (t) we have 


z^{t;s^.) = G(p) 


Re(xi(t)e^*^) 


= G(p)[x^i(t) cos xl(t)sin c/!] (30) 


' cJ^inC^m. 
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since s.,s , and u are assumed independent and have zero mean, 

* M 

E[z^(t;s^.)z^(t;s^O 

K K 

= P2[G(p)6*(p)H-(p)H'*(p)Ry{0)]sin(2c/>) ( 

E[z.^(t;s^)z^{t;Sq)|«?] 

K K 

= -"V PiCG(p)G*(p)H'(p)H'*(p)R^^(0)]sin(2<)!>) ( 

E[z^(t;u)z^{t;u)i‘f] = 0 ( 

since u-j(t), U 2 (t) were assumed independent. Therefore, 

K K 

S{cf) = -1^ G(p)G^p)H'(p)H'*(p)(P 2 Ry( 0 ) - P^R^ (0))sin(2y) 

( 

5.5 Spectral Density of the Equivalent Noise 

Now we will obtain the spectral density of z for fixed ^ . 

We introduce the notation 

Kt) = G(p)H‘ (p)a(t) 

for any signal a(t). Then 

R-{t) = G(p)G*(p)H'(p}H'*(p)R^(x) 

Let C = K, K . Then we have 
1 m 

z^{t;s.) = /C'/Pj sin(«/)y(t) 
z^(t;Sq) = -/C7p:f cos{c/)d^{t) 

Z|.{t;u) = /C^(cos(«/5ui(t)-sin(<;tf)u2(t)) 

Z 3 (t;s.) = /?/Pj cos{<?)y(t) 



J!lnC(. 


om 


where x!(t) and3tl(t) are real. Since 


Similarly, 


x^I{t) = P2H'{p)y(t) 


2 ^(t;s.) = /19^/I^[G(p)H‘(p)y(t}]sin^ 


Z 5 (t;s^.) = -6(p) v^K-j K^' lm( 




= - G(p)[xi(t)sin<P+3$l(t)cos(P] 

" >^'^'^'PJtS(p)H'(p)y(t)]cosc/> 


For the s^{t) effect, we use the fact that 


X (t) = H’ (p)d^(t) 


We find that 


Zc(t;Sq) =’-/iqir/P^CG{p)H''(p)d^.{t)]cos<f (35) 

= /K^y^[G(p)H'(p)d^(t)]sin tP (36) 

For the u(t) effect we have 

z^(t;u) = /Xj 1 TG(p){Re[H'(p)u^(t)+jH'(p)u 2 {t))e'^‘^]} 

= /i^G(p)[H' (p)u-|(t)cos«P -H'(p)u 2 (t)sin(fl (37) 

Zs(t;u) = VK-|K^G(p)[H' (p)u^(t)sin<F+H'(p)u 2 (t)cos ({>] (38) 

Now we can calculate S((p). 


eE[z(t)|<f] = E[z^(t)2^(t) |t^] 


= E[Zj,(t;s.)z^(t;s.)|c^] + E[z^(t;Sq)z^(t;s^) |</>] 
+E[z^(t;u)z^(t;u)ic^] 
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Z5(t;Sq) = /t/p^sin{c^)d^ (t) 


Hence, 


z^(t;u) = /C'(-sin(<j5)u^ (t)-cos(cf>)u 2 (t)) 

S(cp) = E[z(t)|f] = I {P2-Ry(0)-P^R^^(0))sin(2<?) 

1 {z(t) - E[z(t)|^^]) 

= i P2{y^(t)~R-(0))sin(2cf»)- i P^(d2{t)-R^^(0))sin(2^?>) 

- >^P|^ cos(2<f)y(t)d^(t) 

+ /Pj cos(2</>}y{t)u-| (t) - sin( 2 «;<>)y(t}ir 2 (t) 


+ /P^' sin(2^)d^(t)u^ (t) + /Fj' cos{2c(>)d^ {t)‘u’ 2 {t) 

- 1 sin(2c|)u^(t) - cos( 2 r/i)Lr^(t)U 2 {t) 

so +1’ sin(2^0)u2(t} (53 ) 

-4- R^(xl(p). =E{(z(t')-E[z{t) |<f]}(z(t+T)-E[z(t+T) |cf])} 

2 

= \ P2{ECy^(t)y^(t+T)]-R|(0))sin^(2<p)+^^(E[d^{t)d^{t+T)] 

-r|- ( 0))sin^(2ip)) 

+ P^P2 cos^(2<f)R-(T)Rj|T) + P2 cos^(24>)R-(t)R- (t) 

+ P2 sin^(2 (}))R_(T)R- (t) + P^ sin^(2f)Rj^(x)iU^{T) 

+ P^ cos^( 2 (^>)Rj^(t)R- (t) + ^ sin^2^)E[n^(t)u^(t+T)] 

- I sin^(2£())R2 (0) + cos^(2(j>)R^ (t) 

^ n ‘^1 

- 1 sin^(2(^)R|^ (0) + i sin^{2<j!)E[u^(t)u^(t+T)] ( 54) 

This equation for R2 {t|(})) holds because we have assumed LT|{t) and U2 (t) 

are independent and identically distributed for any choice of t and t. 

—2 2 

We must obtain E[u-| (t)u^ (t+x)]. To do this, we define the two random 


cjCinC^ 
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variables X = u-j(t), Y = u^(t+r). Then E[u^(t)u^(t+T)3 = E[X^Y^]. 

E(XYl 22 

Let a = = ayj p = ~ Y • - • To calculate E[X Y ], we introduce the 

a 

zero-mean Gaussian random variable 

W = Y - pX . (55) 

Since E[XW] = E[X(Y-pX)] = 0 = E(X)E{W), then X and W are independent. 
We have 

= E(W^) = E[Y^-2pXY+p^X^] = a^(l-p^) (56) 

Then 

E[xV] = E[X^(W^+2pXW+p^X^)] 

_ 2 2, 2 . 4. . 

- a + p E(X ) 

= /(1-p^) + p^a/ 

= a"^(l+2p2) 

Thus, 

E[Ul(t)Q^(t+T)] = a^1+2p^(r)) 

= ( 0 ) + ZrI ( t ) 

Ui u^ 

so 

" 1 Sin^(2f) [P2Ry2(T)+PiR32(T)] 

0 I 

+ P^P2 cos^( 2 f)R-(T)Rj^(T) + P2Ry(x)R-^{T) 


S.6 Nonlinear Costas Loop Theory 

We are now in a position to investigate the Costas loop performance. 



(57) 


(58) 
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To keep the notation manageable we will develop the theory for a 
first-order loop only. 

The equation of loop operation (12) can be rewritten for a 
first-order loop as 


df(t) 

dt 


% - KyS(if) - 'Kya^{‘^)n^(t) 


(60) 


where n^(t) is a unit-variance zero-mean random process. Approximati 
^ 2 ^^) ^ delta-correlated Gaussian process the above equation can 

be rewritten as an Ito stochastic differential equation 


ng 


d<?(t) = 




1 2 1 
I ^ f tft\ n 


w 


dt 




(61) 


where W(t) is a Brownian motion process. Introducing the notation 

, P cla (f) 

K,W) = Sj - K^,SW) -|p- ( 62 ) 

(drift coefficient) 

K 2 (^) = - Kyo^(c^) (diffusion coefficient) (63) 

the equation takes the form 

dc^(t) = (cf )dt + Kp(f )dW(t) (64) 

and we may use the standard techniques to characterize the 
stationary behavior ofc^. In particular, the p.d.f. of ifi, the 
modulo H tt reduced phase error, is given by 

1 

P(<i>) J* exp UQ(s)ds] (65) 

-TT 

where 


cJ^lnC^c 


om 
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C Is a normalization constant and J is the average rate of 
cycle slips - N_. In the absence of loop stress the density 
function is therefore 

exp£-UQ(<{i)] . (67) 

5.7 Hard-Limiting Shuttle Repeater 

The hard-limiting repeater implementation is illustrated in 
Figure 5.3a. The low-pass filter before the hard-limiter is assumed 
to be wide enough to pass the signal undistorted. 

Assume that signal + noise before the hard-limiter is 

r(t) = i/Fd 2 (t)+ n(t) (68) 

and y(t) is signal + noise after the hard-limiter, i.e., 

y(t) = sgn r(t) (69) 

From Reference 2 we obtain the correlation -function of y(t). 


R^{t). =E[y(t)y(t+T)] 


= E[sgn(r(t)r(t+x))] 
.2 


= erf 


fe) '^ 2 '^ 


) 


+ — exp 

IT 


(^)i: (.)) 


- (-1) (1-Rh (t))] (70) 


where is the kth Hermite polynomial 


c)Cin(^om 


_i di- 
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^2 . n\ V 

— 


^ : 


y 


— 


J 

M(P) i 

LPF 



(a) With Hard-Limiter (b) Without Hard-Limiter 


Figure 5.3. Two Shuttle Repeater Implementations Studied 
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, 2 -u^/2, _ 2 

erf(x) - e ' du - 


-u' 

i ' 

J (\ 


du 


(71) 


Rd (t) = 

°2 


l“hl/T2 » - ^2 

0 i |tI > T, 


(NRZ) 


1 - 3|t|/T 2, M iTg/Z 

M/T2 - 1 , T2/2 < Itj < Tg (BIPHASE) 

0 , |t| 


(72) 


where T£is the pulse duration- of d 2 * If G(p)H’(p) is identity, 

then the equivalent noise correlation function is particularly simple.. 

We have 

R^(tK), = C^[PiP2 cos2(2f)Ry(T)R^^(T) + P2Ry(x)R^^{T) 

\ 

The equivalent noise spectral dens-ity evaluated at zero frequency 
is then 

S2(0|<f) - C^[P,P2 C0S^(2?) J Ry(t)Rj^(t)dt + PjJ R^(t)R^j^(t)dt 

(74) 


Rd/t)R„/t)dt -pj R„^(t)dt] 


where C=K-|K|^,Ry(t) is given in (70), and 




1 

'1 

hi iV^ 

]ji - 1 

T,/2 i hi 

0 , 

hi 


(75) 
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since channel 1 always uses the biphase data format. We use the 
estimate that S 2 (fM = s^{0|cp) for small f. 

5.8 Linear Shuttle Repeater 

Now consider the case where the hard-1 imi ter i s not used and 
M(p) is general 5 illustrated in Figure 5.3b. Then before the 
LPF we have 


r(t) = /P'd 5 (t) + n{t) 


After it we have 


(76) 


y(t) = 


1 




- b(t) + 


1 




v(t) 


(77) 


where 


b(t) = VPM(p)d2(t) 
v(t) = M(p)n(t) 


(78) 

(79) 


and we have scaled y(t) so that E[y (t)] =1. We see that 


Then 


\ir) = PM(p)M*(p)R^^(t) 

\{t) = M(p)M*(p)R^(t) 

" Rjj(0)+Ry(0) 

-2 


For the statistics of y we find 


= r,(0)!r;( 0 T " 2b{t)7(t) . v^t)] 


(80 ) 

(81 ) 

(821 

(83) 


(84) 


' cJCinC^m 
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Ry 2(0 = E[y^t)y 2 (t+T)] - r 2 (o) 


^b( 0 )+Rv( 


lE[b^(t)b^{t+T)] + 2Rg(0)R-(0) 




+ 4Rg{T)R-(T) + E[v^(t)v (t+x)] - (Rb(0)+R-(0))‘^j 

^ 2 [E[b^(t)b^t+x)] + 2Rr{0)R-(0) 

(Rb( 0 )+R/ 0 )r 

+ 4Rb(.T)R-(x) + R^(0) + 2 r|(t) - (Rb(0)+R^(0))^} 

^ 2 fE[b^(t)'5^(t+x)] - r|{ 0) + 4Rb(x)R4x) 

(Rb( 0 )+Rv( 0 )) 

+ 2R?(x)} fe5) 


To obtain E[b^(t)b^(t+x)] - R^(0), we note that 


b(t) = c^q(t-iT) 


where c^- = /P^ or -/F each with probability j and q{t) is the response of the 

G(p)'«'(p)M(p) filter to a data pulse (NRZ or bi-phase) of duration TsT2 
and absolute height 1. 

b^(t) = cX|^q(t-iT)q(t-kT) 

i k 

= p 2D 2] c.C|^q(t-iT)q{t-kT) 

i i kT^i (87 j 


Rg2(x) = ECb2(t)b^(t+T)] - Rg(0) 

= P^\2D q^(t+x-al)2 


c^Cnq(t-iT)q(t-kT)q(t+i 


- JlT) 


i kj^i 


i x\f% 


•q(t+T-nT) P^^E 


(jCinC^k 


fjCinC^c 
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= P^{Z^q^t-iT) X) q^(t+T-£T^ - q^(t-iT)^') 

i i ' i ' 

+ 2P^ (ZE q{t-iT)q(t-kT)q(t+T-iT)q(t+T-kTj^ 

' i k/i 

= P^{zi q^(t-iT) 2] q^(t+T-£T^ - P^^S q^(t-iT^) 
i i i 

+ q{t-iT)q(t+T-iT) ^^-2P^{ S q^(t-iT)q^(t+x-iT^ 

( 88 ) 

where <^»^denotes the average over any time interval of length T-. A 
similar expression holds for R (x), with q(t) replaced by the response 

of the 6(p)H‘(p) filter to a biphase pulse of duration T-|. 


S^(0|<f) = c2jisin2(2<p) 

+ P^ [* R_^(t)dt 

•'—CO 


[R_2(t)+4R^(t)R-(t)+2R^(t)]dt 

b 


"oww (Rb(t)+R7(t))Rj (t)dt 

^ ^ *' ^cn 




Rjj(0)+R- (0) (R|^(t)+Ry(t))R^^(t)dt 

-CO 

•'—CO —CO ' 


where R|j(t) and Ry(t) are given in (80) and (81 ) and C = K^K^- We 
use the estimate that S^(f|f) = S^(0j<p) for small f. To evaluate 
(89) we may use the fact that 


cjCinC^i 
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(90) 


/ Rw^t)R-(t)dt = f |6(f)H'(f)i^^(f)S^(f)df 

•^-co J^co 

where w and x are signals and and are the 'Fourier transforms of 
and R^, respectively. 

5-9 $-Curve Amplitude Plots 

Figures 5.4.through 5,9 are plots of the amplitiTde of the S-curve. 


% = KJrPr CS(<?)/sin(2cp)] 
1 m 2 


G(p)G*(p)H'(p)H'*(p) (r^{ 0) - pi Rj (0)^ 


(91) 


It is plotted as a' function of Shuttle repeater input Ej^/NQ,for channel 2 
having NRZ or biphase pulses and rates of of 2000,192,16Kbps. The channel-! 
data rate is 192 Kbps. Recall that the S-curve is independent of 
the ground station receiver noise. At the top of each figure are 
curves for the case where the repeater is a hard limiter and at 
the bottom are curves for the case of a lowpass-filter repeater. 

Each curve corresponds to a different bandwidth of the combined 
LPF G(f)H*(f) in units of R 2 . 

Following are the assumptions made in obtaining the curves. 

The combined filter G(f)H‘(f) is a one-pole Butterworth filter. 

In the case of HL repeater, the spectrum of this repeater input 
noise is rectangular with a bandwidth of 1.5 R 2 . In the case of 
LPF repeater, the input white noise and signal are both filtered by 
a four-pole Butterworth filter of bandwidth .1 . 5 R 2 . In both cases, 

P-,/P2 = .25. 
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Figure 5.8. for R 2 = 192 Kbps, Channel 2 Biphase, 
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5,10 Variance of Loop Phase Error 

From Reference 3 we have for loop phase error 






(S'(0))‘ 


where Ng/2 = S^(f=0|«p=0) and is the loop bandwidth, when the 
bandwidth of the equivalent noise is much wider than B|^. It can be 
shown that R 2 S^( 0 j 0 ) does not depend on R 2 but only on R^/R 2 » if the 
bandwidth of the G(f)H'(f) filter is taken as a multiple of R 2 * We 
plot 
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in Figures 5.10 through 5.15, using NRZ or biphase pulses and' three 
different data rates for channel 2, Each figure shov/s results for 
both repeaters considered. Assumptions made include those noted 
in the previous section. Additionally, Ej^/Ng = 10 dB for s^ and 
u processes and the spectrum of u is assumed flat over its bandwidth. 

5.11 Filter Bandwidth Selection 

As shown by equation (43), the Costas loop S-curve contains two 
terms whose relative magnitudes depend on the UQPSK power split, the 
payload interrogator bandwidth, the data rates and the IF and arm 
filter bandwidths in the ground station. In the desired mode of operation the 
first term, which depends on the channel 2 statistic, is larger than the second 
term which is determined by the low-rate channel 1, However, if 
the combined IF and arm filter are too narrow they will reduce 
the channel 2 power to the point where the s-econd term is larger 
than the first. In this instance the Costas loop will track the low- 
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Figure 5.10 
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Figure 5,13. ^ for = R2 = 2000 Kbps, Channel 2 Biphase ,E j^/Nq= 10 dB 
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power channel. This problem is particularly serious if the payload 
signal drops out and the payload interrogator receives noise only. 
Assuming a rectangular spectrum for the noise process n (Fig. 5.2) 
the output spectrum is still almost rectangular (Fig. 5.16) with 
the same bandwidth. Since the combined IF and arm-filters do not 
affect the low-rate channel 1 power we can expect the Costas loop 
lock-point to shift by 90 degrees when the arm-filter bandwidth 
is less than approximately 1/4 at the payload interrogator bandwidth 
for a UQPSK power split of 4:1. In the transition region, i.e. 
when the arm filter bandwidth is approximately one quarter of the 
interrogator bandwidth, the loop will either track unsatisfactorily 
or not at all. These conclusions agree with the findings in [4]. 

5.12 Conclusions 

From the plots. Fig. 5.4 to 5.15 it is clear that the payload 
interrogator implementation could be slightly improved, as far as 
subcarrier recovery is concerned, by remvoing the hard-limiter. This 
is borne out by the fact that the resulting S-curve amplitudes are 
higher and less sensitive to the carrier-to-noise ratio at the 
payload interrogator. This small improvement comes at the expense 
of reduced power control for channel 2 which will result in higher 
crosstalk into channel 1. 

For the worst-case conditions (lowest data rate on channel 2, 
Ep/No = 10 dB, i.e. no margin) and for a loop bandwidth of 500 Hz 
the expected rms phase error in the detector will be 5 degrees. 

If the possibility of losing the payload signal exists the relation- 
ship between the payload interrogator bandwidth and the Costas loop 
arm filter bandwidth has to be selected carefully to avoid tracking 






Figure 5.16. Noise Output Spectrum of Hard-Limiter, 
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the signal with a 90 degree offset, 

NOTE : 

In mode 1 of Ku-band return link, the subcarrier is a square 
wave and not a sine wave as assumed in this analysis. The spectrum 
of the modulated subcarrier then is like that of the sine-wave 
subcarrier near 8.5 MHz but has additional harmonics. For the 
case of the LPF repeater, even when R 2 is as large as 2 Mbps, there 
is almost no overlapping of various parts of the spectrum because 
the bandwidth of the LPF is 1.5 which is less than 4,25 MHz. 

When the repeater is a HL, there may be more overlapping when R^ 
is 2 Mbps. Assuming no overlapping, which is true for practically 
all values of R 2 , then as long as the loop arm multipliers are 
insensitive to signals above 4.25 MHz, the arm signals z^ and z^ 
will be the same for square- wave subcarrier as for sine-wave 
subcarrier if the average power of both channels is multiplied 
by //8. 



- 163 - 



jCinCc 


om 


■ REFERENCES 


1. W. C. Lindsey, "Optimum Performance of Costas Type Receivers," 
Axiomatix Report R7502-1 , Feb. 18, 1975. 

2. W. C. Lindsey and W. R. Braun, "Tracking and Data Relay 
Satellite System (TDRSS) Communication Analysis and Modeling 
Study, Phase I Report," LinCom Corp,, prepared under 

NASA Contract No. NAS 5-23591, October 15, 1976, pp. 93-94. 

3. W. R. Braun and W. C. Lindsey, "Carrier Synchronization 
Techniques for Unbalanced QPSK Signals," prepared under 
NASA Contract NAS 5-23591, to be published in IEEE 
Transactions on Communications . September 1978. 

4. Shuttle Ku-Band and S-Band Communications Implementation Study, 
Final Report, Axiomatix Report No. R7902-3, February 28, 1979. 


cJl^inC^vn 




